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Abstract 

A two higgs doublet model with tree level natural flavor conservation can lead to P 
and T odd nucleon-nucleon (NN) interaction through tree level Feynman diagrams. In this 
article we estimate the magnitudes of the parameters that characterize the strength of the 
S-PS type of NN interaction. Stringent constraints on the parameters can be derived from 
the measurement of EDM's of Xe, Cs and Hg atoms in the ground state. We show that 
the predicted values of the parameters lie close to such constraints provided is large. 
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Although in the SM [1] only one higgs doublet was used to give mass to the weak 
gauge bosons and fermions, there is no good theoretical justification for using only one 
doublet. In fact some of the attractive features of the SM like the tree level custodial 
symmetry and natural flavor conservation [NFC] are preserved in suitably chosen versions 
of the two higgs doublet model [THDM]. The minimal supersymmetric standard model [2] 
which has been extremely successful in explaining all the experimental data so far, requires 
two higgs doublets with Y=±l to preserve supersymmetry. 

In this article we shall consider a general THDM in which tree level natural flavor 
conservation is guaranteed by requiring that 0i couples only to RH 7 3 = — | fermions and 
02 couples only to RH 1 3 = | fermions [3] . We shall show that such a model can give rise 
to four fermionic S-PS nucleon-nucleon interaction which are strong enough to produce 
observable P and T odd effects in high precision atomic and molecular experiments. The 
coupling of the neutral higgs bosons 0° and 0° to the Q = | and Q = — \ quarks are given 
by 

L y = -\-& L m lD D^\ + -U^mjuUi^ + h.c.]. (1) 

where we have denoted the Q = 3 and Q = — | quarks generically by U and D respectivvely. 
i and j are the generation indices. Consider first the contribution of Q = | quarks to tthe 
P and T odd NN interaction. We have 

L uu = ^^^(1 - 75)^(1 - 75)^ < > (2a) 

and 

L uu = + 75)^^(1 + 75)^* < 0«* > • (26) 

The P and T odd S-PS mixing term U i UiUj / y 5 Uj will arise from the Im parts of the 
propagators 1 < 0202 > an d * < 02*02* >• Note that % < 0202* > i s rea ^ an d therefore it 
does not contribute to Lpj, v . In unitary gauge [4] we have 0° = | [$1 — ?-^-$3] and 
0° = -jMi- ■ [<E> 2 — i l^M$ 3 ] where $1, $ 2 and $3 are real scalar fields. It then follows that 
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l ptv = - „ 12 1 < $ 2$3 > U i m i uUiU j m jU i'y 5 U j . (4) 



where a summation over generation indices i and j is implied. In order to convert the 
above P and T odd interaction between a pair of U quarks into P and T odd NN interaction 
we have to find the matrix elements of £V UimnjUi and . Ujiiijui^sUj between nucleon 
states [5]. The light quarks u and d, which are the valence quarks for nucleon and whose 
masses are much small compared to the QCD scale A c , contribute very little to the scalar 
matrix element. The heavy quarks s, c, b and t on the other hand exist inside the nucleon 
only as virtual states and therefore their contribution to the scalar and pseudoscalar matrix 
elements can arise only from closed loops. It can be shown that the contribution of each 
heavy quark to the scalar matrix element is given by < N\QhrrihQh\N >= (tunNN 
where Qh stands for a heavy quark and C = 29* ^ ne P seu do-scalar matrix element however 
receives contribution both from light and heavy quark states. The light quark contribution 
to the pseudoscalar matrix element is due to the coupling of mj^u and dtj^d to the 
nucleon via the light intermediate 7To state. They are given by < N\uijsm u u\N >= 
^^m N (-g A )Nij 5 T 3 N and < N\di lb m d d\N >= --JJM— mN (-g A )Ni l5 T 3 N. The 
contribution of each heavy quark to the pseudoscalar matrix element is given by 

< N\Q h i lb m h Q h \N >= ^ mu ~ md m N {-g A )Ni l5 T 3 N - \m N {-g A )Ni lb N. (5) 

2m u +m d 2 

For the Q = | quarks the net contribution to the pseudoscalar matrix element is given by 
[5] 

— 2 777, Tflrl — 

< N\ui~f 5 m u u\N > +2 < N\Q h ij 5 Q h \N > = m N (-g A )Nvy 5 T 3 N 

m u + m d 

-m N (-g A )Ni l5 N. (6a) 



Whereas for the Q = — | quarks the net contribution is given by [5] 

— - 111' — 2/77/^7 — 

< N\di l5 m d d\N > +2 < N\Q h i lb Q h \N > = -m N (-g A )Ni l5 T 3 N 

m u + m d 

-m N (-g° A )Ni l5 N. (66) 

Here g A and g A are the axial vector coupling constants associated with < iV|<9 M J^ 3 |iV > 
and < iV|<9 M J^|iV > respectively. We then get 

< N\L™y\N > = Cm%(-g A ) lVlll ^f 3> NNm l5 N 



~ Cm 2 N (-g A ) 
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NNNvy 5 T 3 N. (7) 



Similarly the P and T odd NN interactions arising from the exchange of neutral higgs 
boson between a pair of D quarks and an U and a D quark are given by 



< N\L™ V \N > = Cm%(-g A ) lV2ll ^ 3> NNm l5 N 
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- CmU-9A) J , ^ ~. | 2 .. NNN1J5T3N. (8) 

and 



< N\L»? V \N > = C<(-^)[ 8< , # f 3> + > ]AfiVjVr ft iy 

\Vi\V \V2\V 

2 i < $ 2 $3 > rn u - 2m d 1 < $i$ 3 > 2m u - m d 

- (m N (-gA) pp + pp ] x 

|^i I v m u + m d \v 2 \v m u + m d 

NNN1J5T3N. (9) 

For q 2 « m\ the propagator factors appearing above can be written as < $i$3 >~ 

~ l Y. n Cl ^ 3n and < *2$3 >~ -*Z) n ■ Here the index n stands for the neutral 

boson mass eigenstates. The orthogonal 3x3 matrix C in connects the gauge eigenstates 
$i with the mass eigenstates $ n . We shall assume that only the lightest neutral higgs (h) 
contributes dominantly to the propgators. In that case we get 
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tNN r Cl h C 3h V C 2 hC3hV , f- 2/ \ AT AT AT AT 

l ptv = —2i — n — 12 H 21 — n — 12 W2G F (m N {-g A )NNNi-f 5 N 

r Cl h C 3h V 1 |f 2 | / Q XX 

~ 7 1 \ 2" ( 1 — r( 2m " ~ + 1 — ^\ m u - 2m d )) 

(m u + m d )m z h \v 2 \ |vi| 2 

(m u + m d )m z h |ui| |v 2 | 2 
NNN1J5T3N. (10) 



Phenomenologically one can write 

^PTV = -C°_ PS ^iViViVz 7 57V - C^^iViViV^raiV. (11) 

where Cly_p S and Cg_ PS are the parameters characterizing the strenghts of P and T odd 
NN interactions [6]. For the two higgs doublet model under consideration they are given 
by 

r»o _ ot( ~Q\ m 2 N v 3 { C lh C 3h C 2h C 3h ] 

C S _ PS - -2C(-^) m 2 ki | k2 |i^-T- + -j^J- ^) 

and 

pi OA / \ m N V [ c lh c 3h 1 |v 2 | , u 

C S _ PS = 2C(-^a) — g - 7 ; r(] — r(2m u - m d ) + , — ^{m u - 2m d )) 

m h \ m u + rn d ) \v 2 \ \vi\ z 

+ < C2hCih x {^—{m u - 2m d ) + T ^(2m u - m d ))]. (13) 
{m u + m d ) \ Vl \ \v 2 \ 2 

Note that the isovector coupling constant Cg_ PS vanishes if we set \vi\ = \v 2 \ and 
m u = m d as required by isospin symmetry. We now turn to the values of the parametrs 
that appear on the rhs of the phenomenological constants Cg_ PS and Cg_ PS . For an 
orthogonal matrix \C\hCzh\ < -5 and \C 2 hC%h\ < -5 . Chiral symmerty breaking exhibited 
in the light pseudoscalar mass spectrum implies that m u = 5 Mev and m d = 7 Mev. 
The direct search for the neutral higgs boson by LEP1 has produced the lower bound 
mh > 58.4 Gev [7] for the three generation SM with minimal Higgs sector and m/ ll > 44 
Gev (m/u < mh 2 ) for the minimal supersymmetric SM. The CDF group associated with 
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the Fermilab tevatron has reported having observed the top quark with a mass of around 
175 Gev [8]. The CDF result implies a large mass heirarchy between the top quark and the 



bottom quark which can naturally arise from a large value of if gt(nit) ~ gb{mb)- In 



the later case we get |fi | ~ 5 Gev and \v 2 \ ~ 175 Gev. The axial vector coupling constants 
g A and g A can be assumed to be approximately equal, i.e. g A ~ g A ~ —1.24. For the above 
values of the parameters we find that \Cg_ PS \ < 3.08 x 10~ 2 and \Cg_ PS \ < 2.30 x 10~ 2 . 

The question that arises is how does the predicted values of the parameters compare 
with the limits derived from experiments searching for P and T violating effects in atomic 
and molecular systems. The P and T odd NN interaction given by eqn. [11] can give rise to 
various T odd nuclear multipoles [6], the simplest of which are the electric dipole moment 
(EDM), Schiff moment (SM) and magnetic quadrupole moment (MQM). The EDM of the 
nucleus does not lead to an overall atomic or molecular EDM because of the well known 
Schiff's theorem [9]. Further for atoms in the ground state with closed electronic shells 
the MQM does not induce an atomic EDM. So the nuclear SM is the only T odd nuclear 
multipole that leads to an atomic EDM. The SM's of 129 < 131 Xe, 133 Cs and 199 ' 201 Hg 
nuclei has been calculated [6] for the P and T odd NN interaction given by eqn. [11]. The 
EDM's of Xe, Cs and Hg atoms have also been related [6] to their respective SM's using the 
relativistic Hartree-Fock technique. The theoretical calculations when combined with the 
measurement of linear Stark effect in the ground state of Xe, Cs and Hg atoms produces 
the following bounds on the isoscalar parameter Cg_ PS [6]: Cg_ PS (Xe) = —.07 ± .24, 
C%_ PS (Cs) = .6 ± 2.4 ± 0.6 and C° s _ PS (Hg) = .01 ± .03. The expeimental limit obtained 
from Hg is the most stringent. The reason being due to larger Z one gains an order of 
magnitude in the EDM by going from Xe or Cs to Hg. Experiments with TIF molecule 
also lead to a bound on Cg_ PS , but it is only an order of magnitude estimate because 
of theoretical uncertainties involved in estimating the SM of Tl nucleus. We find that 
the predicted values of Cg_ PS in our model lie surprisingly close to the experimental 
limits derived from measurements with Hg for which the statistical error is smallest. It is 
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therefore extremely important to improve the accuracy for measuring the EDM of atoms 
and molecules in order to discriminate the different models of T violation on the basis of 
their closeness to the experimental bound and obtain a better physical understanding of 
its origin. 

In conclusion in this article we have estimated the magnitude of the parameters that 
characterize the strength of P and T odd NN interaction in the context of a THDM. The 
model satisfies the conditions of tree level NFC. We find that if the ratio of vev's 



V 2 



then the predicted values of the parameters are close to the best experimental limits ob- 
tained from measurements with Hg. A large value of tan (3 = might naturally explain 
the large t-b mass splitting. However the same could also have other phenomenological 
effects which should be carefully analysed before accepting such a scenario. Present exper- 
iments however do not rule out this possibility and that makes the values of Cg_ PS and 
Cg_ PS presented in this article somewhat interesting. 
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